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Abstract 

Once density-dependent meson masses are introduced into the nu- 
clear many-body problem, conventional mechanisms for saturation no 
longer operate. We suggest that a loop correction, essentially the 
introduction of the axial vector coupling gAiPi k) as function of den- 
sity p and momentum k, can bring about saturation, and present 
schematic calculations to illustrate this. We find that a very small 
density-dependence in gj^ gives rise to a very large saturating effect 
on nuclear matter. In fact, this new saturation mechanism turns out 
to be more powerful than any of the conventional mechanisms. 

PACS numbers: 21.65.+f, 21.60.Jz, 24.85.+p 



For many years during which nonrelativistic Brueckner theory was em- 
ployed, saturation could not be achieved at nuclear-matter density in the 
nuclear many-body problem |l| . The Walecka mean- field theory 0| appeared 
to solve this problem. In this theory, the scalar density p^, which coupled 
to the attractive scalar interactions, increased less rapidly with density than 
the vector density py, which coupled to the repulsive vector interactions. As 
a result, saturation occurred, and this could easily be arranged to happen at 
nuclear matter density po- 

Recent calculations [^ ^, Q have shown that rather complete relativistic 
Brueckner-Hartree-Fock calculations reproduce very much the same scalar 
and vector mean fields from the Bonn potential for the nucleon-nucleon in- 
teraction P, 0] as have been used in Walecka's mean field theory. Thus, the 
mean fields of the Walecka theory can be derived from an interaction which 
describes nucleon-nucleon scattering well. This would seem to bring the story 
of nuclear saturation to a successful conclusion. 

There are, however, new elements in the problem. In a number of papers 
(see, e. g., Ref. Q) it was proposed that meson masses decrease, with in- 
creasing density, at about the same rate as the nucleon effective mass. (The 
quasiparticle velocity vqp = p/m*j^ defines the effective mass.) Reasons for 



the dropping meson masses were most economically summarized by Brown 
and Rho [^ . Rather complete calculations of the density dependence of the 
vector-meson masses m* and m* have now been carried out by Hatsuda et 



al. [|TOl within the QCD sum rule (or finite energy sum rule) formalism. The 



authors find 

^-1-C^ (1) 

my(0) po ^ ' 

where my is the vector- meson mass, with C ~ 0.18. We believe the error 
quoted as ±40% in Ref. [rU| to be larger than the real uncertainty in the 



calculation, which can be cast in a form to depend chiefly on the value of the 
nucleon sigma term Tj^^n'-, which has fairly reliably been pinned down to be 



45 MeV 111 



We wish to point out that one of the most often employed nucleon effective 
masses in nuclear structure calculations is that of the SkM* potential [Q 

(2) 



my 1 



tun 1 + 0.25p/po 
which gives m^^/mN = 0.8 for p = po- Whereas one can see from the 



calculations |jTO[ that the mass of the vector meson rriy mainly follows the 



cube root of the quark condensate < qq >, roughly as 

*\ 1/3 



^-^^i^l , (3) 

mv \< qQ > / 



the nucleoli effective mass involves not only the scalar condensate < qq >*, 
but also the vector one < q^q >, and is much more difficult to calculate in 



the QCD sum rule formalism dTSf. Thus, the m]^/m^ may not accurately 
follow that of my/mv, but it seems reasonable that these ratios are about 
the same at nuclear matter density po! and we shall make the asumption that 
they scale together as in Ref. [||. 



The pion mass m^ does not change appreciably with density fl^. Thus, 



when a common scaling is assumed for meson masses, we must go back and 
make a correction for the lack of scaling of the pion mass. 
Given that masses decrease with a common scale 

A(P) ^ !!^ = !!^ = ... , (4) 

then this common scale can be taken out [§] 

H{r) = \{p)H,ac{\r) = XH^^dx) (5) 

where x = Xr and where H^ac is the original Hamiltonian for the many-body 
system with r replaced by Xr, but with vacuum masses for both nucleon and 
mesons. From Eq. (5), which we assume to hold at least approximately, we 
learn two important results: 



1. There is no saturation, with the assumption Eq. (4), at this level. (Al- 
though fairly obvious, this will be demonstrated below by explicit cal- 
culations.) 

2. Given the scaling Eq. (4), the net binding energy of nuclear matter 
depends only linearly on the effective mass, proportional to m*j^/mj\/, 
a very weak dependence, slightly favoring larger m%f for greater bind- 
ing energies. But we will scarcely be able to calculate the energy E 
accurately enough in order to determine m^ from such a calculation. 

Certainly, the lack of saturation in our theory with scaling masses is a 
grave deficiency. Brown et al. |p forced saturation in an ad hoc fashion by 
including a two-loop correction in the phenomenological theory, but this is 
unsatisfactory in the long run. Since this work, it has been realized |^ that 
the scaling operates at best only at the mean field level, mirroring changes in 
the vacuum condensate. Loop corrections go beyond the scaling Eq. (4) and 
introduce new scales, at least in certain ranges of density. The most familiar 
loop correction in nuclear physics is that which changes the ratio Qa/ Qv from 
unity at the level of tree diagrams (in, e. g., the linear a-model) to 1.26 at 
zero density. It can be calculated by the Adler-Weisberger relation. We shall 



show now that the density dependence of this loop correction can provide a 
mechanism for saturation. 

Whereas (^^ = 1-26 at zero density, there are many empirical indications 



that it decreases with density. Both Rho ||T5[ and Ohta and Wahamatsu [116 
find 
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with b ^ 0.8{gQ)NA, where fnNA is the pion coupling between nucleon and 
A {fnNA ~ 2), and ur is the energy difference between A and nucleon 
{cur ^ 300 MeV). 

There is currently controversy as to whether qaIp) drops this rapidly 
with increasing density, but it is clear that in many experiments Gamow- 
Teller strength is missing [0, and this indicates that qa in nuclei is less 
than gA{0)- 

As noted in the last section, if we begin from the pion coupling at tree 
approximation in the linear a model (which corresponds to the mean field 
approximation of Brown and Rho M), then 



dnNN — {gTvNN)tree 9 A ■ 



(8) 



In the linear a model, the scalar a meson is united, by chiral invariance, with 
the pion in a four-vector. Although, the scalar particle with nis ~ 550 MeV 
used in boson-exchange models is fictions, summarizing some enhancement 
in the correlated two-pion exchange, we shall none the less find it useful to 
treat it symmetrically with the pion, for reasons we explain below. We thus 
give this fictions scalar particle also a loop correction, 

QSNN = {gSNN)tree 9 A ■ (9) 

The precise rate predicted for the decrease, with density, for gA^p) in 
Eq. (6) depends on the value of {g'oJNA, which seems to be in the region []TB 



0.4 < {g'o)NA ^ 0.5 . (10) 

Whatever the precise value of (5'o)7va, gA is predicted to drop from 1.26 
to close to unity as p increases from zero to po- This would mean a large (gi- 
gantic) decrease in scalar strength in the nucleon-nucleon interaction, which 
would obviously upset whatever present agreement there is between calcula- 
tions and nature. 

In order to incorporate the loop correction without upsetting completely 
previous calculations, we must generalize g^ to be a function of both density 



and moment urn 

9a ^ — >gA{p,k). (11) 

As we next show, there is considerable cancellation between the density and 
momentum dependencies, but enough dependence on density will be left for 
it to be important to include qa- 

The two-pion-exchange contribution in the nucleon-nucleon interaction 
gives the lowest-mass effective scalar exchange. Let us consider the density 
dependence in this to lowest order in p; see Fig. 1. Such medium corrections, 
as well as the vr — vr rescattering have recently been calculated in considerable 



detail by Mull, Wambach, and Speth [|T^. Our simple calculations here 
should be viewed as a schematization of the complete calculation, with which 
they agree not only qualitatively, but semiquantitatively. 

The formula Eq. (6) is obtained by retaining only the process of Fig. lb 
for momentum k = 0, but summing the bubble contribution to all orders. 
The renormalized operator gA{p) is then to be used in the nucleon space, 
and the process of Fig. la will be included by configuration mixing. We find 
it more convenient to include both processes Figs, la and lb, in principle 
summed to all orders, in our definition of gA{p, k). 



It is clear that the moinentum dependence of qa cancels most of the 
density dependence, as we now show. The modification in the pion propa- 
gator due to the isobar-hole insertion can, to linear order in the density, be 



expressed through the factor [gO| 



5gA{p, A;) = 1 + 0.8 -^ 
Po 



k"^ + m. 
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where we have left out effects of retardation in the pion propagator. Refer- 
ence [0 takes the typical pion loop momentum to be 



< P >^ 3ml , (13) 

so that the first factor < k'^/{k'^ + ml) >^ 3/4 would seem to predominate 
over (5'o)7VA, giving an increase in qa^P, k). More detailed calculations ||T^, pi 
which include the momentum dependence of {gQ)NA show the latter quantity 
to grow sufficiently with momentum so as to be comparable or slightly greater 
than 3/4 at the momentum given by Eq. (13). It is clear that there is a close 
cancellation between the k"^ /{k"^ + m^) and {gQ)NA terms in Eq. (12). 

In the case of the nucleon particle-hole insertion, the (fi'o(A;))AfAr is found 
to be ~ 15% larger than the < k'^/{k'^ + m^) >^ 3/4 term. (Both are 
weighted by a common Lindhard function in this case.) The net result is 

9 



that qa^P, k), as we define it, is decreased from (7^(0, 0) by a factor such as 
is shown in Eq. (7), but with b replaced by 

beff ^ ^b (14) 

In other words, the near cancellation between repulsive and attractive effects, 
the latter from the momentum dependence, cancel most of the decrease with 
density. In this way, we can understand why investigators have not found it 
necessary to include the loop correction qa to date. The fact that there is 
this large cancellation, leaving a small repulsive effect for the scalar strength, 
is born out by the detailed calculation of Mull et al. W^ . 



In our nuclear matter calculations, we start from the Bonn A potential [0 
and employ the relativistic density-dependent Brueckner-Hartree-Fock for- 
mahsm |^, 0. 

In the bottom of Fig. 2, we show the energy of nuclear matter as function 
of kp for two cases. In Case 1, the nucleon effective mass m^ and meson 
masses mj^ of the Bonn potential, other than rrin (which is held fixed) are 
taken to scale as in Ref. ^ (Table 1, last column, therein). In Case 2 the 
nucleon effective mass is taken to scale as in Eq. (2), and the meson masses, 
other than that of the pion, scale in the same way. The behavior of the 
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curves can be understood from our discussion following Eq. (5). 

The upper curves in Fig. 2, which exhibit saturation, are obtained by 
multiplying gsNN and g^NN by 

^''^'" = 1 + J„P/P0 ''" 

In Case 1, 6e// = 0.035 and in Case 2, befj = 0.045. These are about an 
order of magnitude less than the b ^ {q'qJna of Eq. (7), for the reasons noted. 
Obviously, it will be difficult to calculate b^jf with any precision. 

By comparing the upper curves in Fig. 2 with the lower ones, it is clearly 
seen that even this very small density-dependence of qa has a very large 
saturating effect on nuclear matter. In fact, this effect is even larger than the 
Pauli effect that is known to be the most powerful one among conventional 
saturation mechanisms. 

The saturation energy in both cases shown in Fig. 2 is -15.3 MeV per 
nucleon. The saturation density corresponds to kp = 1.35 fm~^ for Case 1 
and kp = 1-30 fm~^ for Case 2. These saturation densities are better than 
those from relativistic Brueckner-Hartree-Fock calculations where kp = 1.40 
fm~^ is obtained together with an energy of -15.6 MeV per nucleon B. 



In the QCD sum rule calculations [^, ^, ^ the lowest order contri- 
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bution, which gives the main part of the scalar and vector fields, is linear 
in density, just as are the mean fields in the Walecka theory. In this sense, 
the Walecka mean field theory rather directly mirrors QCD p^. Our work 



supports this. The way this comes about in detail, beginning from two-body 
interactions which reproduce the nuclear scattering, such as the Bonn po- 
tential H, 0, is, however, intricate. In Refs. ||, |[ ^ it can be seen that the 
effective scalar and vector mean fields arising from the Bonn potential evalu- 
ated in the relativistic Brueckner-Hartree-Fock (RBHF) calculation are equal 
to the Walecka mean fields only in the vicinity of nuclear matter density. The 
effective coupling constants in RBHF are density dependent, after the effects 
of short-range correlations and exchange are taken into account. For exam- 
ple, at kp = 0.8,1.1,1.5 fm-i Ref. g obtains gl/Ajr = 12.3,8.91,6.23 and 
g^/4:7r = 18.63,13.48,9.06, respectively, in agreement with Ref. p|. At nor- 
mal nuclear matter density, kp = 1.35 fm~^, the scalar field is S = —355.7 
MeV, the vector field V = +275 MeV 0, very close to the values normally 
employed in the Walecka theory. 

For completeness we should note that Drukarev and Levin [^ have cal- 
culated the density dependence oi qa] i- e., gA^p), from QCD sum rules. They 



essentially confirm the results of Refs. |T^, |T6|. We suggest that it may be 

12 



necessary to include the effects of gA{p-, k) in tlie QCD sum rule calculations 
of the binding energy in order to achieve saturation. 

One of the authors (G.E.B.) would like to thank J. Speth and J. Wambach 
for extremely useful discussions. 

References 



[1] 

[2] 
[3] 
[4] 



[5] 
[6] 
[7 



B. D. Day, Comments Nucl. Part. Phys. 11, 115 (1983). 

B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

R. Brockmann and R. Machleidt, Phys. Rev. C 42, 1965 (1990). 

G. Q. Li, R. Machleidt, and R. Brockmann, Phys. Rev. C 45, 2782 
(1992). 

R. Brockmann and H. Toki, Phys. Rev. Lett 68, 340 (1992). 

R. Machleidt, K. Hohnde, and C. Elster, Phys. Rep. 149, 1 (1987). 

R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

G. E. Brown, H. Miither, and M. Prakash, Nucl. Phys. A506, 565 
(1990). 

13 



[9] G. E. Brown and M. Rho, Phys Rev. Lett. 66, 2720 (1991). 

[10] T. Hatsuda and S. H. Lee, Phys. Rev. C 46, R34 (1992). 

[11] J. Gasser, H. Leutwyler, and M. E. Sainio, Phys. Lett. B253, 252 (1991). 

[12] J. Bartel et al, NucL Phys. A386, 79 (1982). 

[13] R. J. Furnstahl, D. K. Griegel, and T. D. Cohen, Preprint OSUPP-92- 
0401. 

[14] G. E. Brown, M. Rho, and V. Koch, Nucl. Phys. A535, 701 (1991). 

[15] M. Rho, Nucl. Phys. A231, 493 (1974). 

[16] K. Ohta and M. Wahamatsu, NucL Phys. A234, 445 (1975). 

[17] See, e. g., C. A. Goodman, "Three Decades of Missing Gamow- Teller 
Strength", in "Spin and Isospin in Nuclear Interactions", S. W. Wissink, 
C. D. Goodman, and G. E. Walker, eds. (Plenum Press, New York and 
London, 1986) p. 167. 

[18] M. Johnson, AIP Gonf. Proc. 163, 352 (1988). 

[19] V. Mull, J. Wambach, and J. Speth, Preprint KFA-IKP(TH)-1992-5. 

14 



[20] T. Ainsworth, G. E. Brown, M. Prakash, and W. Weise, Phys. Lett. 
B200, 413 (1988). 

[21] G. E. Brown, Z. B. Li, and K. F. Liu, Nucl. Phys. A, to be published. 

[22] E. G. Drukarev and E. M. Levin, Nucl. Phys. A511, 679 (1990). 

[23] E. G. Drukarev and M. G. Ryskin, St. Petersburg Nuclear Physics In- 
stitute Preprint (1992). 

[24] T. D. Cohen, R. J. Furnstahl, and D. K. Griegel, Phys. Rev. Lett. 67, 
961 (1991). 

[25] E. G. Drukarev and E. M. Levin, Nucl. Phys. A532, 695 (1991). 



15 



Figure Captions 

Figure 1. Corrections, to linear order, for the density in two-pion ex- 
change. In a) the nucleon particle-hole insertion is shown; in b) the isobar- 
hole insertion. In both cases, a linear density- dependence results from the 
sum over hole states. The cross-hatched insertions in the bubbles represent 
the vertex corrections {g'oJNN and {g'oJNA, respectively. 

Figure 2. Energy per nucleon in nuclear matter. The bottom lines 
(dashed) represent two energies, as function of kp, with scaling masses. 
Curve 1 refers to Case 1 and curve 2 to Case 2 (cf. text). The correc- 
tion < Qa^p) > gives the solid curves in the top of the figure. The dotted hne 
is the conventional relativistic Brueckner-Hartree-Fock (RBHF) result [^]. 
(In all cases, then Bonn A potential is applied.) The cross denotes the 
empirical value for nuclear matter saturation. 
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The figures can be obtained upon request from R. Maclileidt (macli- 
leid@iduil.bitnet; maclileid@tamaluit.pliys.uidalio.edu). Please specify if you 
want them by ordinary mail or telefax and include your mailing address or 
fax number with your request. 
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